Phase-field approach to heterogeneous nucleation 



(N 
O 

o 

(N 

o 

Q 

o 



Mario Castro 
Grupo Interdisciplinar de Sistemas Complejos ( GISC) & 

Escuela Tecnica Superior de Ingenieria (ICAI), 
Universidad Pontificia Comillas, E-28008 Madrid, Spain 

We consider the problem of heterogeneous nucleation and growth. The system is described by 
a phase field model in which the temperature is included through thermal noise. We show that 
this phase field approach is suitable to describe homogeneous as well as heterogeneous nucleation 
starting from several general hypotheses. Thus we can investigate the influence of grain boundaries, 
localized impurities, or any general kind of imperfections in a systematic way. We also put forward 
the applicability of our model to study other physical situations such as island formation, amorphous 
crystallization, or recrystallization. 
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I. INTRODUCTION 

Processes driven by nucleation and growth have at- 
tracted much attention during past decades, from a fun- 
damental point of view [1] and for tailoring some techno- 
logical applications. Some of them are: the recrystalliza- 
tion of deformed metals [2], controlling the nucleation 
and growth of islands on terraces in order to get large 
scale arrays of nanostructures [3], or the manufacturing 
of thin-film transistors which are the basic devices for 
some applications as solar cells [4], random access static 
memories [5], or active matrix-addressed flat-panel dis- 
plays [6]. 

In all the above-mentioned processes, a metastable 
phase decays into a stable one via a fluctuation which 
produces a critical cluster of atoms (for instance, a crit- 
ical island in the case of terrace growth, or a critical 
atom cluster in the case of crystallization). This transi- 
tion is called nucleation. At a certain fixed temperature, 
clusters with sizes greater than a critical one become sta- 
ble nuclei; otherwise they shrink and eventually vanish. 
Such a critical size arises from the competition between 
the surface tension and the chemical potential difference 
between phases, yielding an energy barrier that has to be 
overcome to build up a critical nucleus. For the examples 
presented above, the system can be considered, under 
certain conditions, two dimensional so it is straightfor- 
ward to write the free energy of circular grain of radius 
r: 



AF(r) = 2nra - nr 2 Afj,/n, 



(1) 



where, a, A/i, and Q are the surface tension, the chem- 
ical potential difference between phases, and the mean 
volume occupied by an atom, respectively. 

Notwithstanding, in some practical situations this 
transformation is not perfectly homogeneous due to, for 
instance, the presence of physical boundaries, such as ter- 
race steps, the interplay between different kinds of parti- 
cles, or the appearance of impurities [7], or even to some 
preexisting order embedded in the initial phase formed 



during its manufacturing [8] . The lack of uniformity not 
only catalyzes the transformation but affects the final size 
distribution of grains [9,10]. 

Much work has been devoted to the study of homo- 
geneous nucleation in different contexts [1,11-13], and 
to the so-called site saturation in which nucleation takes 
place just at the beginning of the transformation [14], 
but only a few studies were devoted to the intermediate 
situations in which nucleation is heterogeneous, both in 
space and in time. The spatial extent of heterogeneity 
can be measured using the fraction of transformed mate- 
rial (or more generically, the volume fraction of the newly 
transformed phase) onto the stable phase which often 
obeys the Kolmogorov-Johnson-Mehl-Avrami (KJMA) 
equation [15]: 



X(t) = l- e xp{-[(t-rO/r c ] m } ! 



(2) 



where Ti is the incubation time, t c is the characteristic 
transformation time, and m is an exponent which charac- 
terizes the degree of heterogeneity of the system and its 
dimensionality [9,10]. For instance, for a two-dimensional 
system in the two limiting cases, homogeneous nucleation 
and site saturation, m takes the values 3 and 2, respec- 
tively. 

There have been some approaches to the problem of 
heterogeneous nucleation in the past. For instance, Kar- 
pov et al. [16] generalized the homogeneous case by 
adding random contributions to the surface tension and 
the chemical potential, which affects the nucleation rate. 
A similar static approach was followed by Liu [17], who 
calculated the variation suffered by the surface tension 
and the chemical potential due to the presence of a cir- 
cular impurity. Other authors such as Enomoto [18] 
or Weinberg [19] considered a phenomenological time- 
dependent nucleation rate to analyze heterogeneous nu- 
cleation. More recently, Castro and co-workers [9,10] 
have introduced a lattice model to determine both time- 
dependent and spatial effects of heterogeneities. 

The main aim of this paper is analyze, both analyti- 
cally and numerically, the origin and effect of impurities 
or defects on grain nucleation and subsequent growth. 
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The rest of the paper is organized as follows. In Sec. 
II we introduce a phase field model obtained from a 
functional free energy and report some results for ho- 
mogeneous nuclcation. In Sec. Ill, we include impurities 
through the boundary conditions of the main equation 
and calculate the effect of those impurities in nucleation 
rate and KJMA exponent. In Sec. IV, we consider a more 
general case in which the disorder can be due to mech- 
anisms which are different with the ones presented in 
Sec. III. These mechanisms are introduced as quenched 
noise. Finally, we end the paper by summarizing the 
main results in Sec. V, focusing on the applicability of the 
equations to different physical processes, and discussing 
further generalizations of the model. 



II. EVOLUTION EQUATION 

Phase-field models have been widely studied in the last 
few years, as an efficient computational tool to simulate 
some moving boundary problems which, in the so-called 
sharp interface limit (or sometimes thin interface limit, 
see below), are physically equivalent [14,20-26]. Among 
them, Jou and Lusk [14] studied homogeneous nucleation 
and site saturation using a one-field phase field model. 
The main objection of their approach is the fact that the 
critical clusters are created ad hoc so, on the one hand, 
the model explains the KJMA equation just by construc- 
tion and, on the other hand, it cannot explain the exis- 
tence of an incubation time observed in the experiments. 
A similar approach was used by Roy et al. [12] to study 
nucleation in a phase field model with nonlocal interac- 
tions. Recently, more complex phase field models have 
been proposed for similar systems: The so-called multi- 
phase-field models [27,28], in which every cluster appear- 
ing embedded in the metastable phase is described by its 
own field and, on the other hand, those models in which 
the phase field is coupled with another field representing 
the orientation of crystalline planes [29,30]. The main 
problem concerning the multiphase-field models is that 
the free energy depends explicitly on the grain orienta- 
tion (or phase). This is solved by the second ones but, 
besides this, what depends on the grain orientation is 
the grain-boundary velocity. Although this second kind 
of approach seems to be very promising, in both cases the 
symmetry under rotations of the grain crystalline planes 
is broken. 

In this paper we are interested only in the overall dy- 
namics of the nucleation process, thus, we simply gen- 
eralize the model presented in Ref. [14], supplementing 
it with thermal noise to make explicit the temperature 
dependence of the system. To make clearer this gener- 
alization, we advance that the noise term is the driving 
force for nucleation, so it is not needed to create artifi- 
cially critical clusters every integration time step as in 
Ref.' [14]. 

Let us introduce the main ingredients of the model. We 



define an order parameter 4>, which takes the value — 1 
in the metastable phase and +1 in the stable phase. The 
grain boundary (which separates both phases) is located 
at 4> = [31]. We also define a free-energy functional 
which takes into account the grain-boundary energy and 
the chemical potential difference between phases. Gener- 
ically, we can define 

//IT 2 \ 
dxdz( — |V0| 2 + /(0)-A 3 (0)j, (3) 

where /(0) and (/(</>) are generic functions of the order 
parameter: /(</>) is an even function of (f> with local min- 
ima at ±1 and g(4>) breaks the symmetry between phases. 
As we will see below, W is a typical length scale related 
to the surface tension, and A is a dimensionless parame- 
ter proportional to the chemical potential difference be- 
tween phases. We will assume that the system relaxes 
towards equilibrium according to the following evolution 
equations [32]: 

Td t = W 2 V 2 <P- U^ + Xg^ + e, (4) 

where r is the typical time scale at which the atoms from 
a phase incorporate to the other; and g^ denote the 
partial derivatives of / and g with respect to (p. Finally, 
(9(r, t) is a Gaussian white noise which stands for the ther- 
mal fluctuations of the system, with zero mean and corre- 
lations given by the fluctuation-dissipation theorem [32]: 

(0(r, t)6(r', t')) = 2rk B T5(v - r')5(t - t'), (5) 

T being the temperature at which the transformation 
takes place. 

To be more specific, we choose 

/(0) = -0 2 /2 + ^/4, (6) 

The main advantage of this choice is that Eq. (4) admits 
a simple stationary solution given by: 

^) = -tanh(-^), (7) 

which represents a front of characteristic width W placed 
at z = 0. In the same way, we choose [33]: 

<?(<«=</>- y- (8) 

The main reason to use Eq. (8) instead of the traditional 
one, g((f>) — <fr, is that, in the first case, T has local min- 
ima at ±1 independently of the value of A, otherwise, 
those minima would be A dependent. 

Using Eq. (7) we can make some considerations about 
the stability of a given fluctuation. To compare with the 
classical nucleation theory [34] , let us consider the free en- 
ergy difference between a system which is initially at the 
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mctastable phase, and a circular grain of radius R given 
approximately by fa (r — R). Thus it can be straightfor- 
wardly shown that, in the thin interface limit [22]: 

AT = T[(f> = fa] - T[(f> = -1] ~ 2irRai - nR 2 A^/Q, 

(9) 

where oi — 2^/2W/3 is the surface tension, and A/z/fi = 
X5g, with Sg = g(+l) — g(— 1)- This equation shows the 
competition between the gain arising from the reduction 
of the grain perimeter, and that related to the increasing 
of its size. As we mentioned above, the critical radius 
arises from this competition. Thus if we take into ac- 
count that exp[— AT(R)/kBT] can be interpreted as the 
barrier that has to be overcome to create a nucleus of 
radius R, then the critical radius is the one which mini- 
mizes that barrier, the one which maximizes AT. Hence 



R* = 



2V2W 
3XSg : 



(10) 



and its corresponding critical free energy is given by 

8irW 2 



AT* = 



9X5g 



(11) 



of the free-energy separating this well with the one at 
<f> = +1. Similar results were obtained by Elder et al. [36] 
in the context of eutectic growth. Despite the good agree- 
ment between simulations and theory, it is important to 
stress that the derivation of KJMA equation makes uses 
of some assumptions: infinite system size, uniform nu- 
cleation, spherical particles, and constant growth rate. 
Some care must be taken in this respect. Despite Eq. 
(4) provides circular grains, those grains do not grow at 
constant velocity. Actually, the grain radius is related to 
time through the following equation: 



1 

t= v 



R-R + R* In 



R-R* 
Ro - R* 



(12) 



where Ro is the initial grain radius and V is the grain 
growth velocity which, in the thin interface limit, is given 

by 



V = 



3WX8g 
VSt 



(13) 




FIG. 1. Transformed fraction X(t) vs time for a 500 x 500 
system, (o) numerical integration of Eq. (4) with W = 1, 
r = 0.15, A = 0.6, Ax = 1, At = 0.01, and k B T = 0.1, aver- 
aged over 100 runs. Solid line, KJMA fit with n = 0.75±0.01 
and r c = 1.68 ± 0.01. 

Once we have established the connection between the 
physical system and the model, let us numerically check 
its capability to reproduce the KJMA equation (2). In 
Fig. 1 we show the results from simulation using the Eu- 
lcr integration scheme [35] , and the corresponding KJMA 
fit to the fraction of transformed material, X(t), taking 
m = 3 (measured as the fraction of sites where <f> > 0). 
We have reproduced this result for a wide range of param- 
eters. The first important result concluded from these 
simulations is related to the nonzero value of the incu- 
bation time Tj. This time is almost always present in 
the experiments, and is related to the free-energy dif- 
ference between the well (j) = — 1 and the maximum 



Thus the velocity is constant just after a short time of 
order R* /V . Hence there are slight deviations of the 
numerical data from the KJMA equation up to times 
of order R* /V. Besides this, the infinite size condition 
has not to be fulfilled necessarily. Actually, the KJMA 
equation is still valid when the number of grains con- 
tributing to growth is large, when the system size L, is 
much larger than the characteristic length scale related 
to nuclcation and growth, I = (V/N) 1 / 3 , N being the 
nucleation rate [13]. Finally, the KJMA equation does 
not take into account interfacial effects that govern the 
growth just before the grains meet each other, so there 
are also some differences between the numerical data and 
the KJMA equation predictions, in this case, at the later 
stages of the transformation. Fortunately, the mentioned 
deviations from the KJMA equation at short and long 
times are quite small. Notwithstanding, it is quite con- 
venient to perform the fitting of X(t) between the 1% 
and the 99% to improve it (see Ref. [13] for further de- 
tails about the validity of KJMA equation) . 

To point out the crucial role played by thermal noise 
in Eq. (4), we show in Fig. 2 the spontaneous and con- 
tinuous nuclcation and growth of grains, which are al- 
most circular despite the underlying integration lattice is 
square. In other words, the model captures all the es- 
sential ingredients of a first-order transition, not only in 
terms of the free-energy difference between phases, but 
also in terms of the dynamical path followed from the 
metastable phase to the stable one. 
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FIG. 2. Numerical integration of Eq. (4) with W = 1, 
t = 0.15, A = 0.6, Ax = 1, At = 0.01, and fc s T = 0.1. 
From left to right, from top to bottom, corresponding times 
1.0, 1.5, 2.0, 2.5, 3.0, and 3.5. The darker sites stand for 
<f> ~ — 1 and the brighter ones for <j> ~ +1. It can be clearly 
seen how the initial metastable phase evolves into almost cir- 
cular grains which nucleate and coalesce continuously. 



III. NUCLEATION AT DEFECTS 



In the previous section we have demonstrated the va- 
lidity of our model to describe homogeneous nucleation. 
This section deals with nucleation in the neighborhood 
of some parts of the system which are structurally dif- 
ferent, namely, at domain walls or at the surrounding of 
foreign particles. Thus we will consider two kinds of de- 
fects: Walls and circular defects. Our main assumption 
here is that the particles in the metastable phase do not 
interact with the defects, which is included in the model 
through the boundary condition 




FIG. 3. Numerical integration of Eq. (4) using boundary 
condition (14), with W = 1, r = 0.1, Ax = 1, At = 0.005, 
and ksT = 0.15. The top three figures correspond to A = 0.35 
at times 1, 5, and 10. The bottom three figures correspond 
to A = 0.6 at times 0.50, 0.75, and 1.00. 

As we have mentioned, the other interesting geome- 
try is circular one. Thus Fig. 4 shows how nucleation is 
enhanced at the boundary of a circular impurity. More- 
over, we can use this result to relate the KJMA exponent, 
to, and the effect of impurities using a finite concentra- 
tion c, of circular impurities of small radius. The results 
are plotted in Fig. 5, where time has been rescaled to 
reveal the differences between several to values, for dif- 
ferent temperatures. At this point, we want to stress 
that this dependence of to on temperature has been ob- 
tained qualitatively in some Pdi_ K Sia; crystallization ex- 
periments [38]. Furthermore, if we change the concen- 
tration of impurities, we can also modify the value of to 
(see inset in Fig. 5). This result agrees with that in Ref. 
[10]. Thus the effect of temperature in KJMA equation 
is not only present in the characteristic times, r, and r c , 
but even in the exponent to in a nontrivial way. 



^=V^n| fc = 0, (14) 

where the subscript b stand for boundary; and n is the 
normal coordinate to the defect boundary. 

In order to clarify the relevance of this boundary con- 
dition on nucleation we have integrated Eq. (4) with the 
prescribed condition (14). Figure 3 shows how nucleation 
is 

enhanced at the walls. Actually, this is the most rele- 
vant mechanism of nucleation when the chemical poten- 
tial difference between phases, Afi cx A, is small. The 
situation changes dramatically if we increase A (or if we 
raise the temperature). This can be better understood 
in the context of island formation in which A can be un- 
derstood as the flux of particles arriving at the surface. 
In such case, when the flux of particles is large, islands 
nucleate everywhere in the sample due to the large prob- 
ability of dimer formation [3], as can be seen in Fig. 3. 




1 31) WU 1 su 'M i $ ion 



FIG. 4. Numerical integration of Eq. (4) using boundary 
condition (14) over the circle perimeter, with W = 1, r = 0.1, 
A = 0.5, Ax = 1, At = 0.005, and k B T = 0.1. From left to 
right and from top to bottom, corresponding to equally spaced 
times from 1 to 6. 
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where A measures the degree of heterogeneity of the sam- 
ple. The main part of this generalization is that it does 
not make any specific assumption about the physical sys- 
tem under consideration. Unfortunately, numerical simu- 
lations suggest that it cannot provide a KJMA exponent 
m, different from 3 (homogeneous nucleation) [37]. Nev- 
ertheless, before exploring new possibilities, we will show 
that the quenched white noise given by Eq. (16) unveils 
that disorder enhances nucleation (Fig. 6). 



FIG. 5. Fraction of transformed phase X(t) vs scaled time 
for a 500 x 500 system, with W = 1, t = 0.25, A = 0.8, 
Ax = 1, At = 0.005, averaged over 100 ensemble averages us- 
ing a concentration of impurities c = 0.005, and k B T = 0.01 
(dashed line), k B T = 0.027 (dotted line), and k B T = 0.05 
(solid line). Inset: KJMA exponent m for different values 
of the concentration of impurities c, with W = 1, r = 0.25, 
A = 0.8, Aa; = 1, At = 0.005, and k B T = 0.01. 



IV. GENERALIZED HETEROGENEOUS 
NUCLEATION 
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The results reported in the last section demonstrates 
that our phase field model is a powerful tool to further 
advance in the knowledge and modeling of nucleation and 
growth phenomena, in the homogeneous case and in the 
situation where impurities catalyze the transformation in 
a subtle way. Notwithstanding, the heterogeneities of a 
sample are not always due to isolated impurities. This is 
the main reason for the generalization that we introduce 
in this section. 

We will assume that nucleation is heterogeneous, not in 
a phcnomenological way as in other proposed models [39] , 
but sticking to the classic ideas due to Cahn [40]. Thus 
the system contains regions with some extra energy (for 
instance due to some order produced during deposition of 
the amorphous material), or at which nucleation is more 
probable. Let us show how we can cast this model on a 
mathematical footing. 



A. Quenched white noise 

An intuitive way to introduce the latter idea in our 
model is sustained on the assumption that the chemical 
potential depends locally on position, 



A^ A + p(r), 



(15) 



p(r) being a Gaussian random variable with zero mean 
and correlations given by 



(p(r)p(r'))=A^(r-r'), 



(16) 



FIG. 6. Comparison between homogeneous nucleation 
(A = 0, top panels) and heterogeneous nucleation with 
quenched white noise (A = 0.3, bottom panels), at different 
times (from left to right: t = 20, 50, and 100) with: W = 1, 
t = 1, dx = 1, dt = 0.005, k B T = 0.1, and A = 0.4. 

The fact that white noise is not capable to provide re- 
alistic values of m stems with intuition but, actually, it 
is very physical. Let us make some calculations to show 
it clearly. If we perform again the thin interface limit 
of Eq. (4) with A — > A + p(r), we obtain the stochastic 
version of Eq. (9). Hence: 



AT = 2nRai - nR 2 Ap/n + X (R), 
where x(R) is given by 

X (R) = Sg J dvp(v)e(r - R), 



(17) 



(18) 



being the Heaviside step function. The new noise term 
x(R) has also zero mean and correlations given by: 



( X (R)X(R')) =8gkirR 2 5(R-R'). 



(19) 



For brevity, we denote F = AT, F = 2itR<ti -nR 2 Ap/n 
and D(R) = SgkirR 2 . Consequently, for every value of R 
the free energy is a random number with a distribution 
given by (recall that p is Gaussian) 



P(F) = 



V^D(R) 



-(F-F ) 2 /2D(R) 



(20) 
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Taking into account that exp[— AJ 7 (R)/ k B T] can be un- 
derstood as the barrier that has to be overcome to create 
a nucleus of radius R, the critical radius will be simply 
the one which minimizes that barrier. In the simplest 
case where A = 0, this minimum is the same as the max- 
imum of A.F, but the situation changes because now AT 
is a random variable. Thus in the physical system, we 
can only obtain information about its mean value: 



{e- F ' kBT ) = j dFe- F ' kBT P(F) = exp 



D(R) 



2{k B TY k B T 
(21) 

To be more specific, the critical radius is the solution of 

d (D{R) F \ _ ^ff*_ 2 ^ W 

dR\2k B T °J R=R , 3Sg(X + A/2fc B T)' 

(22) 



This is the result we were seeking: The critical radius 
diminishes so that the nucleation rate increases and also 
depends on temperature. But, as we anticipated above, 
the only effect of quenched white noise is changing the 
time and length scales of the system, but not the nucle- 
ation conditions necessary to obtain m < 3. Moreover, 
if we re-scale both length and time, then the systems 
with A = or A =/= cannot be distinguished. There 
is another interesting point of this theory which should 
be remarked. In some experiments of recrystallization 
the a priori value of AJ-"* is actually much smaller than 
that measured from the experiments [2] , confirming that 
there is some kind of disorder which affects significantly 
this magnitude. 



the white noise case, and also introduced the required 
length scale. As we expected, the KJMA exponent m 
depends on the correlation length l c as we show in Fig. 7 
as also the nucleation rate is larger than in the homoge- 
neous case. These results can be straightforwardly under- 
stood in terms of the involved characteristic length scales, 
the nucleation and growth length scale, I = (V/N) 1 ^ 3 , 
and the quenched noise correlation length, l c . Thus if 
l c <C I then the system hardly see the defects, and it 
is almost homogeneous, so m ~ 3. On the contrary, if 
l c » I, as grains nucleate preferentially at the defects, 
and they grow so fast that they impinge to each other at 
a mean distance l Cl so there is not enough time to allow 
other grains to nucleate. Consequently, m ~ 2 because 
every grain that nucleates does it mainly at the begin- 
ning of the transformation. Note that these assumptions 
are valid whenever both / and l c are much smaller than 
the system size L. Notwithstanding, some care must be 
taken when simulating the model using large values of l c 
because the nucleation events are restricted to the sur- 
roundings of the defects so there are not many grains to 
sustain the validity of the KJMA equation. 




B. Quenched colored noise 



At this point, we must collect some of the successful 
ingredients of the theory and try to find out which must 
be the new ones in order to provide a general model of 
heterogeneous nucleation. Then, the main question is: 
What was implicit in the approach to heterogeneous nu- 
cleation in Sec. Ill which is missing in Sec. IV A? The an- 
swer is related to the new length scale introduced in Sec. 
Ill through the concentration of impurities (and which is 
proportional to c -1 / 2 ). Therefore we need both ingredi- 
ents: disorder and a new independent length scale. The 
simplest way to include both is by means of a quenched 
colored noise term: 

( P (r) P (r'))=Ae-^ 2 /^. (23) 

For the sake of completeness we have performed the simu- 
lations below using Ornstein-Uhlenbeck correlations [35] 
instead of Gaussian as in Eq. (23), but the results did 
not change significantly. With this choice we have the 
same free-energy distribution of probability P(F) as in 



FIG. 7. KJMA exponent, m, for different values of the cor- 
relation length, using W = 1, r = 0.15, A = 0.5, Ax = 1, 
At = 0.05, k B T = 0.05 and A = 0.05. (o) and solid line, 
results from simulations. The dashed lines are a guide to the 
eye. 



Therefore the latter result contains relevant physical 
information relative to the basic conditions yielding an 
exponent in < 3. For instance, it is consistent with those 
crystallization experiments of Sii_a;Gea; by Olivares el 
al. [8] or recrystallization of Pdi_ x Si x by Price [38]. In 
the first case that correlation length could be related 
to the concentration of preexisting SiGe crystals created 
during deposition, and in the second case with the for- 
mer polycrystalline structure which has been deformed 
before recrystallization. 



G 



V. DISCUSSION AND CONCLUSIONS 

In this paper we investigated a phase field model mod- 
eling nucleation and growth events, such as crystalliza- 
tion, island formation, or recrystallization. Starting from 
a Langevin equation we have shown how it is a simple 
(and computationally efficient) way to explore homoge- 
neous nucleation, in terms of the well-known description 
provided by the KJMA equation, Eq. (2). Taking this 
equation as a starting point, we have provided several 
approaches to the problem of heterogeneous nucleation. 
The first approach is related to the existence of defects 
which are inert for the main growing material, such as 
wall domains or impurities. In both cases we have shown 
how nucleation is enhanced near those defects. More- 
over, a finite concentration of isolated defects is able to 
modify the value of the KJMA exponent to, which takes 
values lower than m = 3, the known value for homoge- 
neous nucleation. This kind of defects can be suitable 
to study secondary crystallization, island nucleation in 
heteroepitaxial systems (if one of the growing particles 
presents a slow dynamics compared to the others), or 
nucleation at step edges in terrace growth. In particu- 
lar, wall defects can be used to investigate the influence 
of a terrace step on island formation (note that A plays 
the role of the flux of particles arriving at the terrace). 
In fact, the competition between dimer formation and 
nucleation at step edges has been reported in Ref. [3], 
and can be graphically seen in Fig. 3. Generalizing the 
model to higher dimensions, this kind of wall defect could 
simulate the glass substrate on which the amorphous ma- 
terial is grown prior to crystallization. In the context of 
amorphous crystallization, circular defects can also be 
understood as crystals created at the growing stage, be- 
fore crystallization takes place. 

We have also presented a generalized model for het- 
erogeneous nucleation which deals statistically with im- 
perfections or impurities by means of a quenched noise 
term. The main result is related to the fact that this 
quenched noise must be spatially correlated to provide a 
general description in terms of the KJMA exponent. As 
we mentioned above, the correlation length can be related 
to preexisting crystalline regions embedded in the amor- 
phous host or, some deformed polycrystalline structure, 
but it may also be used to quantify the effect on island 
formation of mechanical stresses caused by the growth of 
multilayer devices. 

Due to the general character of the model, it can be 
easily generalized. For instance, a model that takes into 
account the possible time dependence of the disorder 
should include two fields: the phase-field and another 
one related to the motion of impurities. Besides this, 
sometimes the nucleating grains present some kind of 
anisotropy that should be included in the free-energy 
functional, given by Eq. (3). Moreover, the model pre- 
sented is only suitable for isothermal transformations so, 
to study non isothermal effects, we should couple the 



phase field with an evolution equation for the tempera- 
ture, similar to some prescriptions provided in the con- 
text of non-isothermal cutcctic crystallization [41]. Fi- 
nally, there is some work left to do related to the suit- 
able range of parameters which allow comparisons with 
experimental data and also to improve the efficiency of 
the numerical simulations, in particular when using cir- 
cular impurities to simulate heterogeneous nucleation. 
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